Ancilla-assisted quantum process tomography 
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Complete and precise characterization of a quantum dynamical process can be achieved via the 
method of quantum process tomography. Using a source of correlated photons, we have imple- 
mented several methods, each investigating a wide range of processes, e.g., unitary, decohering, 
and polarizing. One of these methods, ancilla-assisted process tomography (AAPT), makes use of 
an additional "ancilla system," and we have theoretically determined the conditions when AAPT is 
possible. Surprisingly, entanglement is not required. We present data obtained using both separable 
and entangled input states. The use of entanglement yields superior results, however. 

PACS numbers: 42.50.-p, 32.80.-t, 85.60.Gz 



Quantum information science exploits quantum me- 
chanics to achieve information processing tasks impossi- 
ble in the classical world. Recent experiments 2] have 
reported the implementation of a wide variety of simple 
quantum information processing tasks. It is important to 
benchmark the performance of experimental systems as 
quantum information processing devices: one promising 
method, proposed in 1997, is quantum process tomogra- 
phy (QPT) 0. Standard QPT (SQPT) involves prepar- 
ing an ensemble of a number of different quantum states, 
subjecting each of them to the (fixed) quantum process to 
be characterized, and performing quantum state tomog- 
raphy on the outputs. An alternative to SQPT, which we 
refer to as ancilla-assisted process tomography (AAPT), 
introduces an extra ancilla qubit, and involves prepara- 
tion and tomography of only a single 2-qubit quantum 
state, rather than four 1-qubit states 4]. As a special 
case, entanglement-assisted process tomography (EAPT) 
describes the situation when the ancilla is initially max- 
imally entangled with the system being characterized. 

To date, SQPT has been realized in liquid nuclear mag- 
netic resonance systems 5] while SQPT and EAPT have 
been demonstrated in optical systems, but only for uni- 
tary transforms 0. Here we describe optical implemen- 
tations of SQPT, EAPT, and non-entangled AAPT for 
a variety of processes, including unitary, decohering, and 
non-trace preserving (e.g., partial polarizing) operations. 
We also report a theoretical result completely character- 
izing the class of states usable for AAPT. An equivalent 
result was independently developed, and was reported 
just prior to our own 7]. 

In SQPT, a quantum system, A, experiences an un- 
known quantum process, £. To determine S we first 
choose some fixed set of states {pj} which form a basis for 
the set of operators acting on the state space of system 
^4, e.g., {pj} = {pH^ PVi pD^ Pr] for a polarization qubit 
(throughout this paper H, V, D, A, R, and L denote Hori- 



zontal, Vertical, Diagonal, Anti-diagonal, Right-circular, 
and Left-circular polarization, respectively). Each state 
pj is then subject to the process f , and quantum state 
tomography |^ |^, is used to experimentally deter- 
mine the output £{pj). £ is fully characterized if we de- 
termine matrices Ej, known as operation elements, such 
that S{p) = ^jP^ji ^P' This representation is known 
as an operator-sum decomposition . 

In AAPT the process £ is characterized by preparing 
a single state, a, and then measuring {£ 01) (a). This 
requires an ancilla system, B, with Hilbert space dimen- 
sion at least as great as that of A. For an appropriate 
initial state, it is possible to characterize £ by prepar- 
ing the state a, performing the process £ on system A 
— leaving system B completely isolated — and taking a 
tomography of the output {£ 0X){a). The total number 
of measurements is the same in AAPT (16 measurements 
on a single 2-qubit state) as in SQPT (four measurements 
on each of four input states). 

AAPT has advantages over SQPT, most notably being 
that preparation of only a single quantum state is neces- 
sary for its operation. Consider the possibility of using 
it as a diagnostic tool in a quantum computer. When 
an unknown effect acts on less than half of a system of 
qubit s, knowledge of the larger state before and after 
the change is sufficient to exactly predict the effect this 
change will have on every other state. (Assuming that 
the larger state is usable for AAPT — see below). Al- 
ternatively, SQPT has the advantage that it is generally 
easier to produce and measure states with fewer qubits. 

We have investigated a variety of dynamical processes, 
using the three methods of SQPT, EAPT, and non- 
entangled AAPT. Our processes operate on the polariza- 
tion state of a single photon. We used spontaneous para- 
metric downconversion (of a 351-nm pump beam) in a 
nonlinear crystal (BBO) to create pairs of time-correlated 
photons at 702 nm. For SQPT, by triggering on one pho- 
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FIG. 1: Experimental arrangements to perform quantum pro- 
cess tomography. A 351-nm pump is directed through two 
0.6mm-thick BBO crystals, giving rise to pairs of correlated 
photons at 702 nm, which are detected using Si avalanche pho- 
todiodes and fast coincidence electronics. A, B, and C above 
denote which elements are present for SQPT, EAPT, and non- 
entangled A APT, respectively, a.) Single-qubit process to- 
mography: Polarizer (P), half waveplate (HWP) and quarter 
waveplate (QWP) allow preparation of required pure single- 
photon (conditioned on "trigger" detection) states; identical 
elements allow tomography of the post-process states, b.) 
Entanglement-assisted tomography: The source produces the 
maximally entangled state {\HH) — \ VV)) / \f2. A two-photon 
tomography of the output allows reconstruction of the pro- 
cess, c.) Ancilla- assisted tomography: The source produces 
the Werner state pw ~ ^/ + ||7)(7|)5 where I7) is a max- 
imally entangled state. Although there is no entanglement, 
the correlations in pw allow A APT. 

ton, the other was prepared into a single-photon state 
with H polarization (Fig. 1). Half and quarter wave- 
plates converted the horizontal polarization into an ar- 
bitrary state, thus allowing us to prepare the necessary 
input states pu^ Pv-, Pd-, and pR. The tomography of the 
post-process states was performed by measuring (in coin- 
cidence with the trigger detector) the Stokes parameters 
Si = Ph-Pv. S2 = Pd-Pa. and ^3 = Pr-Pl. and per- 
forming a maximum-likelihood estimation of the density 
matrix (Here Pi denotes a probability: calculated as 
the intensity of a state measured in the i^^ basis divided 
by the total intensity.) Typical measurements yielded a 
maximum of 13,000 photon counts over 30 seconds. 

For our EAPT results, two adjacent BBO crystals were 
used to prepare the maximally entangled state |(/)~) = 
{\HH) - \VV))/^/2) [la]. One of the resulting qubits 
was subjected to the given process, and two-qubit to- 
mography of the pair was then performed by measuring 
the polarization correlations of the photons with 16 mea- 
surements, e.g., in the following bases: HH, HV, HD, 
HR, VH, VV, etc. 0. Note from Fig. 1 that the ele- 
ments used in SQPT to prepare the single photon state 
are now placed (in reverse order) in the other detection 
arm, highlighting the symmetry of the two techniques. 

We also performed AAPT using the non-entangled 
Werner state pw = where I7) is a max- 

imally entangled state. To prepare this state we ad- 
just the polarization of the pump beam until the down 



conversion crystals produce the pure, partially entan- 
gled state 1^^) + i- A half wave- 
plate at 22.5° in each arm then transforms this state 
into 10) = ^\HH) + ^^\HV) + ^^\VH) + ^^\VV). 
Next we pass each photon through a decoherer, an 11-cm 
piece of quartz which separates the horizontal and verti- 
cal components of the polarization by ~ 100 /im, which 
is the coherence length of the individual photons (de- 
termined by the 3-mm diameter collection irises and the 
5-nm bandwidth (FWHM) interference filters.) This de- 
stroys all coherence terms in except for \HH){yV\ 
and \VV){HH\. An additional, shorter decoherer in the 
idler arm lowers these terms to achieve a state which has 
99.2 ± .8% fidelity 13] with the above Werner state. 

For single-qubit processes, a convenient graphical rep- 
resentation plots the transformation of the sphere of all 
possible states (e.g., the Poincare sphere for polariza- 
tion) m, as determined by the action of the process on 
the set of basis states, pj. For example, all unitary trans- 
formations are equivalent to a rotation about some axis 
(Fig. 2b). Decoherence is represented by a collapsing 
of the sphere toward a "spindle" (Fig. 2c); for instance, 
complete decoherence in the HV basis leaves the states 
\H) and \V) unmodified, but transforms the states \D) 
and \R) into the completely mixed state at the center of 
the sphere. This graphical approach can even be applied 
to lossy processes, e.g., partial polarizers, though it is 
important to note that it does not indicate the amount 
of loss, only the quantum state of the surviving qubits. 

We now outline the general procedure for SQPT, as 
described in 3]. Rather than directly determining the 
operation elements Ej^ SQPT relates these to a fixed 
set of operators, {Em}^ where Ej = ^^e-jmEm and 
ejm can be complex. This allows us to define a sin- 
gle matrix, Xi that fully characterizes thej)rocess: if 
we rewrite the process as £{p) = ^rnn^rnpElxmn then 
X is a positive Hermitian matrix, Xmn = Y^k^kmel^- 
See Fig. 3 for examples of experimentally determined 
X matrices. To determine x from a set of measure- 
ments, we choose a set of basis states {pj}, such that for 
each input state pj, state tomography returns an output, 

^{Pj) = Hk^jkPk' If we define EmPjE\ = T.kf^Tk'Pk 
(where /3Jl^ is another complex number matrix which 
we determine from our choice of input basis states {pj}, 
output basis states {p/c}, and operators {Ej}), we can 
see that Y.kT.mnXmnPjr pk = Y^k^jkPk, independent 
of pk] P is invertible; and Xmn = Y^jk {f^~^)Jk ^jk- 

In our experiment we use {Em} = {I ,(Tx,(Jy,(Jz}i re- 
spectively equivalent to the following optical elements: 
nothing; a HWP plate at 45°; an optically active ele- 
ment; a HWP plate at 0°. The diagonal elements of the 
X-matrix correspond, respectively, to the probability of 
carrying out the I,(Jx,(Jy, and processes, while the 
off-diagonal elements correspond to coherence processes 
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FIG. 2: Geometric mappings for three quantum processes - 
(a) identity , (b) unitary transformation, and (c) decoherence 
- measured using SQPT (left), EAPT (center), and AAPT 
(right). The axes are the Stokes parameters {Si, Ss)- 
The colored mesh surfaces show how all pure states are trans- 
formed by the process. The initial states H, R, V, and A are 
shown by the green, red, yellow, and blue dots, respectively. 
The transformation of initial mixed states (inside the surface) 
may be interpolated from the pure state results using the lin- 
earity of quantum mechanics. The mesh coloring denotes the 
orientation of the transformed sphere. 

of the form axpcy and (jypcjx^ etc. 

We investigated several processes, including the iden- 
tity, a unitary rotation, a decoherer, and both a coher- 
ent and an incoherent partial polarizer (see below). The 
results for the identity process measure how well the in- 
put state(s) are preserved. We used SQPT, EAPT, and 
AAPT to measure the same unitary rotation process (a 
birefringent waveplate). The results were in close agree- 
ment (Fig. 2b); the resulting x matrices had an aver- 
age process fidelity Q between the three methods, of 

= 100.4 ± .8%. Likewise, the SQPT and EAPT mea- 
surements of a decohering process (implemented with a 
6.3-mm piece of quartz) yielded = 99. 9 ±.3% (Fig. 2c). 
The same process, when measured using our Werner 



Real Imaginary ' Real Imaginary 




FIG. 3: x-i^citrices determined from EAPT for (a) unitary 
and (b) decohering processes, as shown in Fig. 2. 



State, appears to be a recoherer — a process which is not 
a positive map. Recall that this Werner state was pre- 
pared using a thick piece of quartz to temporally separate 
the H and V components of the light, introducing deco- 
herence. Consider adding another piece of quartz, with 
optic axis perpendicular to the first, after the original. 
This also temporally shifts the H and V components of 
the light, but in the opposite direction, undoing the orig- 
inal decoherence. Our decohering process does exactly 
this, effectively recohering the Werner State — impossi- 
ble for a 1-qubit process. The resolution to this paradox 
lies in the implicit assumption that the measured process 
does not act on any degrees of freedom used to prepare 
the input state other than the tested qubit. For example, 
if frequency is traced over to prepare a mixed input state, 
a process that couples to frequency cannot be measured. 

Coherent and incoherent partial polarizers were mea- 
sured in order to highlight the role coherence plays in 
lossy processes. A glass plate at Brewster's angle to an 
incident beam is a coherent partial polarizer, as the op- 
eration of the plate maintains the pre-existing phase re- 
lationship between the horizontal component of the light 
(completely transmitted) and the vertical component of 
the light (partially reflected). For the incoherent case, 
consider inserting a horizontal polarizer into the beam 
50% of the time. Half the time only the horizontal com- 
ponent of the light will be transmitted, but more impor- 
tantly, the transmitted light will have no coherence re- 
lationship with the light that does not pass through the 
polarizer. For the coherent partial polarizer, pure states 
remain pure but slide toward H along the surface of the 
sphere. In the incoherent case pure states travel linearly 
through the sphere to H, becoming mixed (Fig. 4). 

What class of initial states a of the AB system may 
be used for AAPT? This question can be answered us- 
ing an operator generalization of the Schmidt decom- 
position for entangled states Q. To explain this de- 
composition, we introduce an inner product on oper- 
ators, (M^N) = tT{M^N)^ and define an orthonor- 
mal operator basis to be a set of operators {Mj} such 
that (Mj.Mk) tr(MjMfe) Sjk. (For example, an 
orthonormal basis for single-qubit operators is the set 
{I/V2,ax/V2,ajV2,az/V2}). The operator- Schmidt 
decomposition [13 states that an operator M acting on 
AB can be decomposed as M = "^iSiAi Bi, where 
the si are non- negative real numbers, and the sets {Ai} 
and {Bi} form orthonormal operator bases for systems A 
and 5, respectively 1^. The Schmidt number Sch(M) of 
an operator M is defined as the number of non-zero 
terms in the Schmidt decomposition. 

A state a of AB may be used to perform AAPT if 
and only if the Schmidt number of a is where dA is 
the dimension of the state space of system A. Consider 
that in order to measure the mapping of the entire space, 
the input state must possess correlations - represented by 
the Schmidt number - between enough states to form a 
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FIG. 4: Geometric mappings and x matrices for (a) coherent 
and (b) incoherent partially polarizing processes. The former 
was implemented using two glass plates (microscope slides) 
near Brewster's angle [Th ^ 88%, Ty 45%]. The latter was 
simulated by inserting a horizontal polarizer 50% of the time. 
(Real components shown; imaginary contributions < 1%.) 

basis for the mapping. To prove this, expand a in its 
Schmidt decomposition as a = ^iSiAi Bi. Assume 
a has Schmidt number so that the Ai form an or- 
thonormal operator basis, and si > for all I. Let a' be 
the output obtained after letting £ act on system A, that 
is, a' = (£01) {a) = siS{Ai) ^ Bi. By the orthonor- 
mality of the Bi and the previous equation it follows that 
tiB{{I^Bl)a') = si£{Ai)tT{BlBi) = Sm£{Am), and 
so £{Am) = trB((/ BI^)(7') / Sm- The fact that the 
Schmidt number of a is d\ ensures that 5m > 0, so there 
is no problem with division by zero. It follows that it is 
possible to determine the action of £ on an operator basis 
by doing state tomography on <j', and applying the above 
equation. The techniques described earlier can then be 
used to generate a x matrix or transformed sphere. 

Conversely, let Ea be the space of trace-preserving 
quantum operations on system ^4, and let Sab be the 
space of quantum states on system AB. Define a map 
f ■■ Ea ^ Sab by f{£) = {£®I){a). For AAPT, 
we require that / be a one-to-one map, i.e., there are 
never two distinct operations such that f{£i) = /(^2)- 
A parameter counting argument shows that / cannot be 
one-to-one when a has Schmidt number less than d\. 
The dimensionality of the manifold Ea is d\ — d\. Since 
f{£) = si£{Ai)^Bi^ the dimension of the image man- 
ifold /{Ea) is at most Sch(M) x (d^ — 1), because the map 
£ £{Ai) has image of dimension at most — 1. Thus, 
for AAPT we require that Sch(M) x {d\ -l)>d\- d% 
which is only possible when Sch(M) = d\. 

Note that AAPT is possible only when the dimension 
of system B is at least as great as the dimension of system 
A. When this is true, almost all states of system AB may 



be used for AAPT, because the set of states with Schmidt 
number less than d\ has measure zero. That is, a maxi- 
mally entangled input is not required for AAPT — indeed 
many of the viable input states are not entangled at all, 
as demonstrated by our Werner state AAPT. However, 
while almost any state can be used for AAPT, maximally 
entangled states appear to be experimentally optimal in 
that they have perfect non-local correlations. Fig. 2 high- 
lights this difference, as the AAPT results have signifi- 
cantly greater statistical errors than the EAPT (both 
were from identical measurement runs). This compara- 
tive usefulness of entangled versus separable states was 
first introduced and is discussed further in 0. 
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